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There is strong experimental evidence for pairing of polaronic carriers in the normal state, two 
distinct energy scales, d-wave superconducting order parameter, and charge segregation in the form 
of stripes in several cuprates. All these remarkable phenomena might be unified in the framework of 
the bipolaron theory as a result of the formation of mobile bipolarons in the normal state and their 
Bose- Einstein condensation. Extending the BCS theory towards an intermediate and strong-coupling 
regime we show that there are two energy scales in this regime, a temperature independent incoherent 
gap and a temperature dependent coherent gap combining into one temperature dependent global 
gap. The temperature dependence of the gap and single particle (Giaver) tunnelling spectra in 
cuprates are quantitatively described. A framework for understanding of two distinct energy scales 
observed in Giaver tunnelling and Andreev reflection experiments is provided. We suggest that both 
d-wave superconducting order parameter and striped charge distribution result from the bipolaron 
(center-of-mass) energy band dispersion rather than from any particular interaction. 
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Introduction 

There is strong evidence for normal state gaps in high-T c 
cuprates from the uniform susceptibility [QJ|], inelastic 
neutron scattering and NMR electron energy-loss 

spectroscopy ||, specific heat H, angle-resolved photoe- 
mission (ARPES) @,||, tunnelling §,0, photoexcited 
quasiparticle relaxation p| , and some kinetic measure- 
ments One view supported by ARPES is that the 
gap reflects precursor superconducting correlations in the 
BCS-like state below some characteristic temperature T* 
without long range phase coherence. Testing of this hy- 
pothesis with specific heat || and tunnelling |)| data 
find that it cannot be sustained. In particular, there 
is no sign that the gap closes at a given temperature 
T*. On the other hand the strong-coupling extension 
of the BCS theory based on the multi-polaron perturba- 
tion technique firmly predicts the transition to a charged 
Bose liquid in the crossover region of the BCS coupling 
constant A ~ 1 (Bi)polaronic nature of carriers 
in cuprates, supported by the infrared spectroscopy and 
quite a few other experiments |14|, i n particular by the 
isotope effect on the carrier mass |15| provides a natural 
microscopic explanation of the normal state gap ||l6|| . In 
the framework of the bipolaron theory the ground state of 
cuprates is a charged Bose-liquid of intersite mobile bipo- 
larons where single polarons exist as excitations with the 
energy A ~ T* or larger |l7| . A characteristic tempera- 
ture T* of the normal phase is a crossover temperature 
of the order of A where the population of the polaronic 
band becomes comparable with the bipolaron density. 
Along this line the normal state kinetics of cuprates has 
been explained and a theory of tunnelling and 
angle-resoleved photoemission [|l9| in cuprates has been 
developed describing essential spectral features of a single 
hole in these doped Mott insulators. 

However the temperature and doping dependence of 
the gap still remains a subject of controversy. More- 
over, reflection experiments, in which an incoming elec- 
tron from the normal side of a normal/superconducting 
contact is reflected as a hole along the same trajectory 
f^of , revealed a much smaller gap edge than the bias at 
the tunnelling conductance maxima in a few underdoped 
cuprates |2l[| . Recent intrinsic tunnelling measurements 
on a series of Bi '2212' single crystals 0] showed dis- 
tinctly different behaviour of the superconducting and 
normal state gaps with the magnetic field. Such coex- 
istance of two distinct gaps in cuprates is not well un- 
derstood ||^,^]. There is also a d-like superconducting 
order parameter (changing sign when the CuOi plane is 
rotated by 7r/2) in cuprates as observed in the phase- 
sensitive experiments [ p4[ . And, finally, there is growing 
experimental evidence |25| , p(| that stripes occur in doped 
cuprates. The C11O2 plane of high-T c cuprates has been 
found to be inhomogeneous on the microscopic scale, sug- 



gesting that holes condense into large stripe domains [^7) . 

Here I suggest that two distinct energy scales, d-wave 
superconducting order parameter, and charge segrega- 
tion in cuprates are manifestations of one and the same 
phenomenon, which is the normal state pairing of pola- 
ronic holes in the doped charge-transfer Mott insulators. 
The main assumption is that a short-range attraction 
potential in cuprates is large compared with the (renor- 
malised) polaron Fermi energy. Our main point is in- 
dependent of the microscopic nature of the attraction. 
Real-space pairs might be lattice and/or spin bipolarons 
p7f ; or any other preformed pairs. 

2. Normal and superconducting gap in cuprates 



Recently we have proposed a toy model p8[ , which ac- 
counts for two energy scales. The model is a one- 
dimensional Hamiltonian including kinetic energy of car- 
riers in an effective mass (m) approximation and a local 
attraction potential, V(x — x') — —US(x — x 1 ) as 

H = J2jdxiPt(x) (-±.±p-^ 

S V / 

-U J dx^ix^lix^^x^ix), (1) 

where s =f , j is the spin (Ti = ks = 1). The first band to 
be doped in cuprates is the oxygen band inside the Hub- 
bard gap as established in polarised photoemission Jl9| . 
This band is almost one dimensional as discussed in Rcf. 
]l8[ , so that our (quasi) one-dimensional approximation 
is a realistic starting point. Solving a two-particle prob- 
lem with the 5-function potential one obtains a bound 
state with the binding energy 2A p = ^mU 2 , and with 
the radius of the bound state r = 2/mU. We assume 
that this radius is less than the inter-carrier distance in 
cuprates, which puts a constraint on the doping level, 
Ep < 2A p , where Ep is the polaron Fermi energy. Then 
real-space pairs are formed. If three-dimensional cor- 
rections to the energy spectrum of pairs are taken into 
account (see, for example, Ref. [|9|) the ground state 
of the system is the Bose-Einstein condensate (BEC). 
The chemical potential is pinned below the band edge by 
about A p both in the superconducting and normal state 
]j~7|] , so that the normal state single-particle gap is A p . 

Two distinct energy scales appear naturally if we take 
into account that in the superconducting state (T < T c ) 
the single-particle excitations interact with the conden- 
sate via the same potential U . Applying the Bogoliubov 
approximation [ |30f we reduce the Hamiltonian, Eq.(l) to 
a quadratic form as 

* = £/^I(*)(-^-m)^) 
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dx[A c ip\(x)^\{x) +H.c 



(2) 



where the coherent pairing potential, 
A c = —U (ipiix^^x)), is proportional to the square root 
of the condensate density, A c = constant x no(T) 1 / 2 . 
The single-particle excitation energy spectrum E(k) is 
found using the Bogoliubov transformation as 



E(k) = [(fc 2 /2m + A p ) 2 + A 2 ] 



1/2 



(3) 



This spectrum is quite different from the BCS quasipar- 
ticles because the chemical potential is negative with re- 
spect to the bottom of the single-particle band, u = — A p . 
The single particle gap, A, defined as the minimum of 
E(k), is given by 



A = [A 2 + A 2 ] 



I 1/2 



(4) 



It varies with 

1 /2 

temperature from A(0) = [A 2 + A c (0) 2 ] at zero tem- 
perature down to the temperature independent A p above 
T c . The condensate density depends on temperature as 
1 - {T/T c ) d / 2 in the ideal three (d = 3) and (quasi) two- 
dimensional (d — 2) Bose-gas. In the three-dimensional 
charged Bose-gas it has an exponential temperature de- 
pendence at low temperatures due to a plasma gap in 
the Bogoliubov excitation spectrum Q , which might be 
highly anisotropic in cuprates ]l7|| . Near T c one can ex- 
pect a power law dependence, no(T) oc 1 — (T/T c ) n with 
n > d/2 because the condensate plasmon |H| depends 
on temperature. As shown in Ref. |^8| the theoretical 
temperature dependence, Eq.(4) describes well the pio- 
neering experimental observation of the anomalous gap 
in YBa2Cu30j^s in the electron-energy-loss spectra by 
Demuth et al f|, with A C (T) 2 = A c (0) 2 x [1 - (T/T c ) n \ 
below T c and zero above T c , and n = 4. 



conductance, Eq.(5) with one of the best STM spectra 
measured in TVi-substituted Bi2Sr2CaCu20g,+ x single 
crystals by Hancottee et al |32|, Fig. la. This experiment 
showed anomalously large 2A/T C > 12. 

The theoretical conductance, Eq.(5) describes well the 
anomalous gap/T c ratio, injection/emission assymmctry, 
zero-bias conductance at zero temperature, and the spec- 
tral shape inside and outside the gap region. There is no 
doubt that the gap, Fig.l is s-like. The conductance, 
Eq.(5) fits also well the conductance curve obtained on 
'pure' Bi2212 single crystals |52|, while a simple d-wave 
BCS denstity of states cannot describe the excess spec- 
tral weight in the peaks and the shape of the conductance 
outside the gap region. We notice that the scattering 
rate, e, derived from the fit, is apparently smaller in the 
'pure' sample than in Ni- substituted one, as expected. 

The same toy model provides a simple theory of tun- 
nelling into bosonic (bipolaronic) superconductor in the 
metallic (no-barrier) regime |ps| . As in the canonical 
BCS approach applied to the Andreev tunnelling by 
Blonder, Tinkham and Klapwijk (3^], the incoming elec- 
tron produces only outgoing particles in the supercon- 
ductor (x > I), allowing for a reflected electron and (An- 
dreev) hole in the normal metal (x < 0). There is also 
a buffer layer of the thickness I at the normal metal- 
superconductor boundary ( x = 0), where the chemical 
potential with respect to the bottom of the conduction 
band changes gradually from a positive large value u in 
the metal to a negative value — A p in the bosonic super- 
conductor. We approximate this buffer layer by a layer 
with a constant chemical potential pi, (— A p < /x& < fi) 
and with the same strength of the pairing potential A c as 
in the bulk superconductor. The Bogoliubov-de Genncs 
equations may be written as usual p3| , with the only 
difference that the chemical potential with respect to the 
bottom of the band is a function of the coordinate x, 



3. Giaver tunnelling and Andreev reflection 

The normal metal-superconductor (SIN) tunnelling con- 
ductance via a dielectric contact, dl/dV is proportional 
to the density of states, p(E) of the spectrum Eq.(3). 
Taking also into account the scattering of single-particle 
excitations by a random potential, thermal lattice and 
spin fluctuations as described in Ref. one finds at 
T = 



dl/dV oc [p 



2eV - 2A 



, ,'-2eV-2A\. 
Ap[ ; ], (5) 



with 

(r ,_ 4 Ai{-2tj)Ai'{-2Z) + Bi{-2Z)Bi\-2Q 
P ^> it 2 [Ai{-2^Y + Bi{-2iYY 

A is the asymmetry coefficient |L8|], Ai(x), Bi{x) the Airy 
functions, and e is the scattering rate. We compare the 



-(l/2m)d 2 /dx 2 
A n 



p(x) A c 
(l/2m)d 2 /dx 2 



p(x) 



tp(x) 
E^{x) (7) 



Thus the two-componet wave function in the normal 
metal is given by 



tp n (x < 0) = | Q | < 



1 > " lq+x -J- b { I | e~ lq+x + a(° l | < ' 



(8) 



while in the buffer layer it has the form 



^,(0 <x<l)=a( e*P +x +p( e-*" x 

V E+S J V E-£ J 

+ 7 ( jL ) e- zp+x +s( 2^_) e^ x , (9) 

V E+£ J \ E-i J 

where the momenta associated with the energy E are 
q ± = [2m(p ± E)] 1 ' 2 and p ± = [2m(p b ± il)} 1/2 with 
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£ = (E 2 — A 2 ) 1 / 2 . The well-behaved solution in the su- 
perconductor with negative chemical potential is given 
by 

j> s (x >l) = c( t.) e lk+x +d( t ) e ik ~ x , (10) 
V £+? / V E-i J 

where the momenta associated with the energy E 
are fc ± = [2m(-A p ± £)] 1/2 - Tne coefficients 
a, b, c, d, a, /3, 7, <5 are determined from the boundary con- 
ditions, which are continuity of ip(x) and its first deriva- 
tive at x = and x = I. Applying the boundary condi- 
tions, and carrying out an algebraic reduction, we find 

a = 2A c q+(p+f-g+ -p-f+g-)/D 1 (11) 

b = -1 + 2q+[(E + 0/ + (<rr - p~g-)/D 

-2q+(E-t)f-(q-f + -p + g + )]/D, (12) 

with 

D = (E + 0(q + f + +P + g + Kq-f- ~ P~ 9~) 

~(E-0(q + .r +p-g-)(q-f + -P + g + ), (13) 

and = p^ cos(p ± l) — ik sin(p ± Z), g ± — k^ cos(p ± l) — 
ip^ s'm^^l). 

The transmisson coefficient for electrical current, 1 + 
\a\ 2 — \b\ 2 is shown in Fig. 2 for different values of I when 
the coherent gap A c is smaller than the pair-breaking 
gap A p . We find two distinct energy scales, one is A c in 
the subgap region due to electron-hole reflection and the 
other one is A, which is the single-particle band edge. 
We notice that the transmission has no subgap structure 
if the buffer layer is absent (I = 0) in both cases. In the 
extreme case of a wide buffer layer, I >> (2mA p ) -1 / 2 , 
Fig. 2. there are some oscillations of the transmission due 
to the bound states inside the buffer layer. It was shown 
in Ref. [j35| that the pair-breaking gap A p is inverse pro- 
portional to the doping level. On the other hand, the 
coherent gap A c scales with the condensate density, and 
therefore with the critical temperature, determined as 
the Bose-Einstein condensation temperature of strongly 
anisotropic 3D bosons [p9| . Therefore we expect that 
A p >> A c in the underdoped cuprates, Fig. 2. Thus the 
model accounts for the two different gaps experimentally 
observed in Giaver tunnelling and electron-hole reflection 
in the underdoped cuprates . 

One of the important experimental findings is the zero- 
bias conductance peak (ZBCP) observed in a number 
of in- plane tunnel junctions |34]. Within our approach, 
ZBCP-like feature appears in the junctions with A p > A c 
and I >> 1 (see thin dashed line in Fig. 2). Overdoped 
cuprates might be in the BCS or intermediate regime 
because the increasing number of holes leads to an in- 
crease of their Fermi energy and, at the same time, to a 
decrease of A p due to screening p5[ . Then, the condi- 
tion for BEC, Ef < 2A P no longer holds, and the d-wave 



BCS approach |}6| might account for ZBCP as well. Such 
evolution from the BEC regime in underdoped cuprates 
to the BCS-like regime in the overdoped samples has 
been found in the conserving T-matrix approximation 
of the two-dimensional attractive Hubbard model fl37| . 
The binding energy of real-space pairs monotonically de- 
creased with increasing density, and vanished at a critical 
density of holes n cr — 0.19 (per cell). 

4. D-wave Bose-Einstein condensate 

To account for the symmetry of the order parameter 
and the orientation dependence of ZBCP our simple one- 
dimensional continuos model can be readily extended to 
two and three dimensions, discrete lattices, and nonlocal 
pair potential. Than, quite generally, the symmetry of 
the single-particle spectrum, Eq.(3) and the symmetry 
of the Bose-Einstein condensate wave-function (i.e. of 
the order parameter) are not necessary the same in the 
bosonic superconductor |Q. The single-particle spec- 
trum is dominated by the incoherent (s-wave) gap, A p , in 
agreement with the c-axis tunnelling data, Fig.l. On the 
other hand, the symmetry of the Bose-Einstein conden- 
sate depends on the bosonic pair (center of mass) energy 
band dispersion. 

Different bipolaron configurations can be found with 
computer simulation techniques based on the minimiza- 
tion of the ground state energy of an ionic lattice with 
two holes. The intersite pairing of the in-plane oxygen 
hole with the apex one is energetically favorable in the 
layered perovskite structures as established by Catlow et 
al J?8|| . This apex or peroxy-like bipolaron can tunnel 
from one cell to another via a direct single polaron tun- 
nelling from one apex oxygen to its apex neighbor. The 
bipolaron band structure has been derived in Ref. |39[ as 

E k V = tcos(k x ^ y ) - t'cos{k ytX ). (14) 

Here the in-plane lattice constant is taken as a = 1, t is 
twice the bipolaron hopping integral between p orbitals 
of the same symmetry elongated in the direction of the 
hopping (ppcr) and —t' is twice the hopping integral in 
the perpendicular direction (ppn). These hopping inte- 
grals are proportional to the single-particle hopping in- 
tegrals between apex oxygen ions. The bipolaron energy 
spectrum in the tight binding approximation consists of 
two bands E x,v formed by the overlap of p x and p y apex 
polaron orbitals, respectively. The energy band minima 
are found at the Brillouin zone boundary, k = (±7r, 0) 
and k = (0, ±7r) rather than at the T point owing to the 
opposite sign of the ppa and ppix hopping integrals. Nei- 
ther band is invariant under crystal symmetry but the 
degenerate doublet is an irreducible representation; un- 
der a 7r/2 rotation the x band transforms into y and vice 
versa. 

If the bipolaron density is low, the bipolaron Hamil- 
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tonian can be mapped onto the charged Bose-gas [ p"7| . 
Charged bosons are condensed below T c into the states 
of the Brillouin zone with the lowest energy, which are 
k = (±7r, 0) and k = (0, ±7r) for the x and y bipolarons, 
respectively. These four states are degenerate, so the 
condensate field-operator \&(m) in the real (site) space 
m = (m x , m y ) is given by 



*±(m) = N~ 1/2 



E 

k=(±7r,0),(0,±7r) 



bkexp(— ik ■ m). (15) 



where N is the number of cells in the crystal, and &k is 
the bipolaron (boson) annihilation operator in k space. 
This is a c- number below T c , so that the condensate wave 
function, which is the superconducting order parameter, 
is given by 



^±(m) = rt^J 2 [cos(7rm x ) ± cos^m^)] , 



(16) 



where n c is the number of bosons per cell in the conden- 
sate. Other combinations of the four degenerate states do 
not respect time-reversal and (or) parity symmetry. The 
two solutions, Eq.(16), are physically identical being re- 
lated by: ^ + (m Xl m y ) — ^ ^(m x ,m y + 1). They have 
(i-wave symmetry changing sign when the CuOi plane is 
rotated by 7r/2 around (0,0) for or around (0,1) for 
The e?-wave symmetry is entirely due to the bipo- 
laron energy dispersion with four minima at k ^ 0. 



5. Bipolaronic stripes 

The antiferromagnetic interaction are thought to give rise 
to charge segregation (stripes) in cuprates j4C],|4l|]. At 
the same time the electron-phonon interaction is strong 
in ionic cuprates and manganites [§~7],(l2) . While for wide 
band polar semiconductors and polymers the charge seg- 
regation (strings) were discussed some time ago |Q, a 
role of the long-range Frohlich electron-phonon interac- 
tion in the charge segregation in narrow band ionic in- 
sulators like cuprates has been addressed only recently 
ffi] , |45f with the opposite conclusions about existance of 
strings. 

We have proved jl5) that the Frohlich interaction com- 
bined with the direct Coulomb repulsion does not lead to 
charge segregation like strings in doped narrow-band in- 
sulators, both in the nonadiabatic and adiabatic regimes. 
The polarisation potential energy of M fermions trapped 
in a string of the length N, 



U 



M 2 I 



N ■ 



(17) 



appears to be lower in value than the repulsive Coulomb 
energy because K — eo£oo/ (eo — £oo) is always larger than 
the high-frequency dielectric constant eoo. Here eo is the 
static dielectric constant, and the integral Ijy is given by 



TV 



(2tt) 3 



7T />7T />7T 

dx dy dz 



sin(iVa;/2) 2 



_ T "J^ iV2 s in(x/2)2 



x (3 — cos x — cos y — cos z) . (18) 
It has the following asymptotics at large N 
1.31 + In N 



In = 



N 



(19) 



The Coulomb energy of spinless one-dimensional 
fermions comprising both Hartree and exchange terms 
is M 



En = 



e 2 M (M - 1) 



[0.916 + In M], 



(20) 



so that the polarisation and Coulomb energy per particle 
becomes 



U/M 



e 2 M 



[0.916 + hiM-a(l. 31 + ln7V)], (21) 



where a = 1 — £ooAo < 1- Minimising this energy with 
respect to the length of the string TV we find 



N = M 1/a exp(-0.31 + 0.916/a), 



(22) 



and 



(U/M) r 



-M x - 1/a exp(0.31 - 0.916/a). (23) 



One can see that the potential energy per particle in- 
creases with the number of particles. Hence, the en- 
ergy of M well separated polarons is lower than the en- 
ergy of polarons trapped in a string no matter corre- 
lated or not. The opposite conclusion derived in Ref. 
p4| originates in an erroneous approximation of the in- 
tegral ijv N 0A5 /N. The correct asymptotic result is 
I N = ln{N)/N. 

In prinicple, the situation might be different if the anti- 
ferromagnetic interaction of doped holes with the parent 
Mott insulator is strong enough |^,^l| . Due to the long- 
range nature of the Coulomb repulsion the length of a 
single stripe should be finite (see, also Ref. |2^,^]). In- 
cluding the Frohlich, Eq.(17) and Coulomb Eq.(20) con- 
tributions one can readily estimate its length for any type 
of the short-range intersite attraction, E att as pq| 



N 



exp 



e aE att Suj 



2.31 



(24) 



where uj is the characteristic frequency of bosonic excita- 
tions (like magnons) responsible for the short-range at- 
traction, Slo its maximum dispersion, and a is the lattice 
constant. Then, with typical values of eo = 30, a = 3.8A 
and E at t = 2 J = 0.3 eV of the t — J model one can hardly 
expect any charge segregation due to antiferromagnetic 
spin fluctuations as well, because N < 2 according to 
Eq.(24). Hence it is more probable that the lattice and 
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spin polarons in cuprates remain in a liquid state in the 
relevant region of the parameters. 

The Frohlich interaction is not the only electron- 
phonon interaction in ionic solids. As discussed in Ref. 
p7[ , any short range electron-phonon interaction, like, for 
example, the deformation potential and/or Jahn- Teller 
distortion can overcome the residual weak repulsion of 
Frohlich polarons to form small bipolarons. At large dis- 
tances small bipolarons weakly repel each other due to 
the long-range Coulomb interaction, strongly suppressed 
by optical phonons. Hence, they form a liquid state in the 
relevant region of their densities and masses p8| , which 
is also confirmed by recent study jl9) of the Holstein- 
Hubbard model. The ground state of the ID Holstein- 
Hubbard model is a liquid of intersite bipolarons with 
a significantly reduced mass |^9|. The bound states of 
three or more polarons are not stable. 

If (bi)polaronic carriers in cuprates are in the liquid 
state, one can pose a key question how stripes can be seen 
at all. Here I suggest that the bipolaron liquid might be 
striped owing to the bipolaron energy band dispersion, 
Eq.(14), rather than to any particular interaction. If 
bipolaron bands have their minima at finite k (includ- 
ing the Brillouin zone boundaries), then the condensate 
wave function is modulated in the real space Eq.(16). As 
a result, the hole density, which is about twice of the con- 
densate density at low temperatures, is striped, with the 
characteristic period of stripes determined by the inverse 
band-minima wave vectors. If disorder is strong enough, 
the condensate might be localised within the wide and 
shallow random potential wells, so bipolaronic stripes can 
also exist in the insulating (due to disorder) oxides. If the 
condensate is striped, one can expect rather different val- 
ues of the coherent gap, A c , for single-particle excitations 
tunnelling along (110) and (100), and, hence the different 
Andreev reflection spectra. Our interpretation of stripes 
is consistent with recent inelastic neutron scattering stud- 
ies of YBa2Cu307_a [|5(J , pj]| , where the incommensurate 
peaks in neutron scattering have been only observed in 
the superconducting state. The vanishing at T c of the 
incommensurate peaks can be actually anticipated over 
a wide part of the phase diagram and in other neutron 
data |p2[. It is, of course, inconsistent with the usual 
stripe picture where a characteristic distance needs to 
be observed in the normal state as well. In our picture 
with the d-wave striped Bose-Einstein condensate the in- 
commensurate neutron peaks disappear at T c or slightly 
above T c (due to superconducting fluctuations), as ob- 
served. 

To conclude I suggest that two distinct energy scales, 
d-wave superconducting order parameter and stripes in 
underdoped cuprates are manifestations of mobile bipo- 
larons. The normal state gap in the charge channel is 
half of the bipolaron binding energy. The superconduct- 
ing gap comprises additional coherent part owing to the 
interaction of single holes with the Bose-Einstein con- 



densate, d-wave superconducting order parameter and 
stripes are intimately connected as the result of the bipo- 
laron energy band dispersion with the minima at finite 
wave vectors of the center-of-mass Brillouin zone. 

I greatly appreciate enlightening discussions with An- 
tonio Bianconi and Robert Laughlin. 
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Figure Captures 

Fig.l. Theoretical tunnelling conductance, Eq.(5) 
(line) compared with the STM conductance in Ni- 
substituted Bi 2 Sr 2 CaCu 2 O s<+x Q (dots) with 2A = 90 
meV, A = 1.05, e = 40 meV. 

Fig. 2. Transmission versus voltage (measured in units 
of Ap/e) for A c = 0.2A p , [i = 10A p , fi b = 2A p and / = 
(thick line), / = 1 (thick dashed line), I — 4 (thin line), 
and I = 8 (thin dashed line) (in units of l/(2mA p ) 1 / 2 ). 
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